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The  elTects  of  ussteai^  sOain  rate  on  the  burning  velocity  of  hydrogen/air  premixed  dainee  have  been 
studied  in  an  opposed  ooale  ooofiguraQoD.  Tbe  numerica]  method  employs  adaptive  time  islegraOoo  of 
a  system  of  differenOal-algebraic  equations.  Detailed  hydrogen/ air  lonetic  mechanism  and  Lraosport  prop* 
erties  are  considered.  The  eqiuvaJeooe  ratio  is  varied  Irom  lean  to  nob  premixtures  in  order  to  change  the 
elective  Leviis  number.  Steai^  Markstein  numbers  for  small  sDain  rate  are  computed  and  compared  with 
experiment  Different  deSnitions  of  dame  burning  velocity  are  examined  under  steady  and  unsteady  flow 
conditions.  It  is  found  that,  as  the  unsteady  Irequency  increases,  large  deplations  between  different  flame 
speeds  are  noted  depending  on  the  location  of  the  flame  speed  evaluation.  Unsteady  flame  response  is 
invesO^ted  m  terms  of  tbe  Markstein  transfer  function,  which  depends  on  the  frequency  of  osciflation 
In  most  cases,  the  flame  speed  variation  attenuates  at  higher  frequencies,  as  tbe  unsteady  frequency 
becomes  comparable  to  the  inverse  of  the  characteristic  flame  time.  Furthermore,  unique  resonance*li]ce 
behavior  is  observed  for  a  range  of  rich  mixture  conditions,  consistent  with  previous  studies  with  linearized 
theory. 


Introduction 

In  tbe  kminar  flamelet  reeime  of  turbulent  pte- 
mixed  combustion,  the  turbulent  burning  velocity  Is 
determined  by  the  integral  of  the  local  laminar 
flamelet  speed  ov^  the  entire  fiame  front  whidi  Is 
corrugated  by  tbe  turbulent  eddies.  It  is  well  known 
that  ^e  bminar  flame  speed  is  not  a  constant,  but  a 
fujiction  of  flame  stretch,  consisting  of  strain  and 
curvature  components  induced  by  the  turbulent 
flow.  There  are  also  flame  instability  modes  that  con¬ 
tribute  to  an  enliajiced  surface  area.  Tlierefbre,  it  is 
importajit  to  understand  ajid  parameterize  the  effect 
of  stretch  on  tbe  flame  speed  in  accurately  predicting 
the  overall  burning  rate  of  premixed  combustion  de¬ 
vices. 

Througli  a  number  of  theoretical  and  experimen¬ 
tal  studies  over  the  past  decades,  it  has  been  fairly 
well  estabhshed  that  tbe  flame  speed,  5,  is  Hnearly 
proportional  to  stretch  in  the  weak  stretch  limit.  This 
relation  is  writtea  In  nondimensional  form  as 

S/Sl  »  1  -  MaKa  (1) 

where  Sx.  Is  the  unstretched  laminar  fiame  speed,  Ka 
is  the  Karlovltz  number  representing  stretch,  and 
Ma  is  tbe  proportionality  constant  known  as  the 
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Markstein  number  [L],  Asymptotic  analysis  showed 
that  the  Markstein  number  is  a  function  of  ther¬ 
modynamic  and  transport  properties  of  the  mixture 
[2^],  such  that  It  can  assume  positive  or  negative 
values  depending  on  tbe  diihisivity  of  tbe  deficient 
reactant  species,^own  as  tbe  Levv^  number  effect. 
Elxperimental  measurements  have  been  performed 
to  confirm  this  prediction  (5-8],  while  computa¬ 
tional  studies  with  detailed  chemical  kinetics  and 
transport  provided  means  to  extract  accurate  values 
of  tbe  Markstein  numbers  for  various  fuel  mixtures 
(9,10). 

Most  of  the  previous  studies  regarding  tbe  stretch 
effect  on  premixed  flames,  however,  are  based  on 
quasi-steady  flames  ajid  bence  may  become  inaccu¬ 
rate  in  hi^y  transient  flows  produced  by  turbu¬ 
lence  .  The  effect  of  unstea<ty  flow  on  premixed  flame 

Firopasation  was  investigated  in  asymptotic  studies 
11,12],  where  it  was  found  tliat  tbe  flame  speed  be¬ 
comes  independent  of  strain  in  the  liigb  frequency 
limit,  wliile  tbe  curvature  component  remains.  Re¬ 
cent  direct  numerical  simulations  of  turbulent  meth¬ 
ane/air  [13]  and  hydrogen/air  [14]  premixed  flames 
also  revealed  that  the  flame  speed  response  becomes 
less  sensitive  to  turbulent  flows  for  higher  turbu¬ 
lence  intensities,  Tbese  results  clearty  suggest  that 
flow  unsteadiness  Is  an  important  parameter  In  mod- 
itying  the  premixed  flame  response. 

In  this  paper,  we  present  a  parametric  study  on 
tbe  response  of  premixed  flames  with  detailed  by- 
drogen/air  chemisby  to  unstea<ty  strain  rate.  Effects 
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of  thermodiifusive  imbaknce  are  investigated  by 
comparing  the  results  for  different  equivalence  ra' 
tios.  Various  definitioas  of  tbe  fiame  speed  arising 
from  finite  fiame  structure  are  also  discussed.  The 
results  are  compared  with  theoretical  predictions 
and  some  recent  experimental  measurements. 


Formuintiun  and  N'umerical  Method 

The  numerical  configuration  Is  a  counterflow  sys¬ 
tem  betsveen  two  opposing  axisymmetrlc  nozzles 
separated  by  a  distance  of  L.  Following  tlie  previous 
work  (IS.lfiX  a  modified  system  of  equations  for  the 
unsteady  problem  is  adopted  In  order  to  take  into 
account  the  high-index  nature  of  the  problem, 
thereby  facilitating  uumerical  Integration  using 
adaptive  time  integration.  Defining  ^e  self-similar 
variables  for  the  velocity,  scaled  radial  velocity, 
temperature,  and  species  mass  fractions  as 

b  »  u{t.  x).  v/t  -  V(f,  x) 


T  »  r(f,  x),  y*  »  Yiit.  x)  (2) 

and  introducing  the  acoustic  pressure,  p  »  P  —  po, 
where  P  is  the  total  static  pressure  and  po  is  the  ther¬ 
modynamic  pressure,  the  conservation  equations  are 
written  as  mass  continuity: 


p  ^ 

P  * 


T  dt  ?  W),  A 


—  (pu)  +  2pV  »  0 

dx 


axial  momentum: 


du  du  dp  dV  4  d  /  du 

0 —  +  />u  —  +  —  —  %i  —  —  — \u  — 

dt  dx  dx  dx  3ai\  dx 


(3) 


radial  momentum; 


(4) 


dV  dV 

p —  +  pH —  + 


dt 


dx 


I 


-—(u'—]  +  A-0 
dx 


(5) 


energy  conservation; 


ar  aT 


species  conservation; 

p  ^  +  pu  ^  ^  -  Wicoj, 

A  tfc  dr 

»  0.  Jc  »  I,  •  -  • ,  K  (7) 

Tlie  equation  of  state  becomes 


p  »  (p  +  Po)W/RT  (8) 

In  the  above  equations,  A{t)  *  (l/r)(0p/ar)is  the 
radial  pressure-gradient  eigenvalue,  p  Is  the  viscos¬ 
ity,  Cp  Is  the  specific  heat  of  the  mixture,  X  is  the 
thermal  conductivity  of  the  mixture,  is  the  en¬ 
thalpy  of  species  x.  Wit  is  the  molecular  weight  of 
species  K,  and  a>i,  Is  the  molar  reaction  rate  for  spe¬ 
cies  K.  The  above  system  of  differential-algebraic 
equations  is  Integrated  using  DASPK  [17],  which 
employs  a  backward-difference  formula  method 
with  adaptive  time-step  and  order  control  The  code 
interacts  with  CHEMKIN  [18]  and  TRANSPORT 

[19]  packages  for  computing  reaction  rates  and 
transport  properties.  The  kinetic  mechanism  used  is 
the  hydrogen/alr  system  developed  by  Yetter  et  al. 

[20] . 

For  tlie  boundaiy  conditions,  fi’esh  reactajits  at 
300  K  are  supplied  at  x  ^  L,  and  burned  product 
was  supplied  at  x  »  0.  where  the  temperature  and 
composition  of  the  product  is  determijied  from  the 
downstream  values  of  the  freely  propagating  pre¬ 
mixed  fiame  for  the  given  fresh  mixture  condition  at 
X  »  L,  In  this  study.  L  Is  set  to  I  cm.  Pressure  Is  set 
to  1  atm  througliout  the  study.  While  the  hydrogen/ 
air  system  behavior  is  known  to  be  significantly  af¬ 
fected  by  pressure,  as  (ar  as  the  effects  of  stretch  on 
the  fiame  speed  are  concerned,  we  expected  that  the 
results  will  be  qualitatively  similar  at  different  pres¬ 
sures.  as  can  be  deduced  from  the  experimental  re¬ 
sults  for  various  pressures  [7].  Tlie  scaled  radial  ve¬ 
locity  V  is  set  to  zero  at  both  boundaries,  and  the 
unsteady  flow  field  Is  imposed  by  specifying  the 
equal  axial  velocity  at  both  boundaries  as  a  (unction 
of  time.  A  sinusoidal  velocity  of  the  form 

u(r  »  0)  »  — u(x  »  L) 

»  «o(l  +  A[1  -  cos{2;^)j}  (9) 

is  imposed  such  that  the  velocity  oscillates  from  uo 
to  Uo  (1  +  2A)  at  a  frequency  of  /  IIz,  which  is  a 
parameter  characterizing  the  unsteady  flow  time 
scale. 

The  spatial  differencing  of  the  equations  follows 
the  previous  study  [16],  The  grid  dependence  of  the 
numerical  solution  was  tested  as  the  grid  refinement 
was  made,  for  both  steady  ajid  unstea<fy  calculations. 


Various  Definitions  of  the  Flame  Speed 

While  the  mathematical  definition  for  the  flame 
speed  ajid  strain  rate  is  based  on  an  asymptotically 
^bin  reaction  zone,  in  reality  the  flame  has  a  finite 
thickness  such  that  these  quantities  cannot  be 
uniquely  defined.  In  majiy  counterflow  experiments 
[21]^  the  fiame  speed  has  l)een  measured  at  the  lead¬ 
ing  edge  of  the  preheat  zone,  which  is  identified  as 
the  local  maximum  of  the  axial  velocity  profile.  This 
upstream  velocity,  denoted  as  herein,  is  relatively 
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Fig  1.  Steady  response  of  varioui 
flaine  speeds  to  oormalized  strain 
rale  Kiior  ^  *  04. 


easy  to  measure  using  laser  Doppler  velndmetiy,  A 
recent  experimental  stu^  of  ujisteady  premixed 
flames  (22]  also  followed  tills  convention.  As  will  be 
discussed  later,  however,  under  unsteady  conditions 
such  a  definition  can  lead  to  significantly  dliforent 
results  compared  to  other  definitions  as  the  unsteady 
frequency  Increases. 

In  numerical  studies  with  finite  rate  chemistry 
more  mathematically  rigorous  definitions  have  been 
suggested  [23].  The  displacement  speed,  derived 
from  the  Hamilton'Jacobi  equation,  is  defined  as  the 
velocity  of  the  flame  surface  represented  hy  an  iso> 
level  surface  of  a  chosen  species.  To  minimise  ther' 
mal  expansion  effects  across  the  flame,  the  density' 
weiglited  displacement  speed,  5d,  for  species  k  Is 
defined  as  [24] 


p,  l?y*l  IVY*I 


3  /  SYi 

_  p/x  _5 


pDk 

Pa 


(^-n) 


(10) 


where  Is  the  density  of  the  upstream  reactajitmlX' 
hire,  Dj.  Is  the  mass  diffosivity  of  spedes  k,  and  tj  Is 
the  coordinate  normal  to  tlie  flame  front.  The  last 
term  in  eq^tion  10  represents  the  efiect  of  flame 
cujvature."V  •  o,  where  n  »  —  Vy*/l7Yj,l  is  the  unit 
vector  normal  to  the  front.  One  disadvantage  of  this 
definition  is  that  the  displacement  speed  depends  on 
the  choice  of  the  isolevel  surface,  although  it  was 
observed  that  the  variation  Is  relatively  small  when 
the  Isolevel  surface  is  near  the  reaction  zone  [13]. 

The  consumption  speed  is  also  defined  [23]  as  a 
measure  of  the  consumption  rate  of  spedes  k  per 
unit  flame  area,  written  as 


"  />u(Ytb  -  VtJ 

where  subscripts  u  ajid  b  denote  the  unbumed  and 
burned  sides,  respecti\«ly.  Tliis  quajitlty  is  direct^ 
related  to  the  overall  burning  rate  and  can  be  of 
practical  importance. 

The  primary  focus  of  this  paper  is  the  effect  of 
strain  rate  on  tlie  flame  speed.  For  the  axisymmetric 
counterflow,  the  Karlovitz  number  represents  the 
strain  component  only,  ajid  is  defined  as 


where 


(12) 


«  =  (13) 

ajid  6  is  the  flame  thickness.  Note  tliat  the  Karlovitz 
number  varies  significantly  in  magnitude  depending 
on  the  definition  of  the  flame  thickness.  We  clioose 
S  B  D1(/Sl  where  Dj,  is  tlie  mass  diffosivity  of  species 
k  with  respect  to  the  upstream  conditions  for  mix¬ 
ture  and  temperature.  For  the  results  presented 
herein,  O2  is  chosen  for  the  marker  spedes. 


Steady  Response 

As  a  reference  point,  we  computed  the  steady  re¬ 
sponse  of  premixed  dames  of  hydrogen/air  mixture 
with  various  equivalence  ratios.  For  each  mixture 
condition,  we  compared  various  definitions  of  the 
flame  speed  described  in  the  previous  section.  As  an 
example.  Fig.  1  shows  the  result  for  a  lean  mixture 
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Fig.  2.  Maj43leis  Dumber  based  on  various  flame  speed 
definitions  as  a  functioo  of  eqiuvaleoce  ratio 


TABLE  1 

Parameters  for  various  test  cases 


Sl 

Do, 

uo 

A 

SJh. 

atuo 

K© 

04 

22.15 

02263 

100 

05 

1.553 

346 

0.5 

59.89 

02317 

140 

0.25 

1.105 

1220 

07 

148.50 

0242L 

500 

0.3 

1006 

1567 

3.0 

258.93 

0.3355 

600 

0.3 

0.957 

1781 

5.0 

130.38 

0.3915 

240 

025 

0.872 

256 

65 

66.40 

04243 

120 

05 

0.315 

505 

3.0  (diluted) 

33.^ 

0  2707 

80 

05 

0.255 

344 

Note.  Units  are  cm/s  for  Sl  and  oo.  cmVs  for  Dog,  and 
s~‘  for  Vq.  In  ad  cases,  the  air  is  uodJuted,  0^{0^  +  N^) 
*  0^1,  eKcept  for  the  last  case  denoted  dduUd  uhere 
(V(Oj  +  Nj)  *  0.114. 

of  s  0.4.  The  solid  lines  denote  the  displacement 
speed  evaluated  at  four  different  values  of  the 
mass  fraction,  where  the  manmum  heat  release  oo 
curs  appnwimatety  at  Yo,  =  0  ,17.  Altliough  quali' 
tabvely  similar,  the  four  values  lead  to  quantitatively 
different  results  for  5^.  The  consumption  speed, 
udiich  is  uniquely  defined  for  a  given  species,  falls 
within  the  S ^  curves  with  the  same  qualitative  trend, 
while  Sy,  the  flame  speed  based  on  the  location  far 
upstream,  tends  to  significantly  overestimate  the  re> 
suit.  The  strain  rate  given  by  equation  13  Is  com* 
puted  at  the  chosen  isocontour  for  Since  no 
spedflc  Isocontour  is  needed  to  compute  and  5^, 
the  strain  rate  value  at  —  0.17  Is  used  here. 

Results  for  the  Marlcsteiu  uumher  obtained  from 
a  linear  curve  fit  for  various  flame  speed  definitions 
are  plotted  in  Fig.  2.  Tlie  range  of  Kj  Is  limited  to 
smaller  values,  approximate^  Ka  <  0,2  for  tlie  def' 
Inition  based  on  Do,.  The  experimental  results  by 


Aung  et  al  [7j  are  also  converted  and  plotted  Quan* 
titative  discrepancies  may  be  attributed  to  the  dif* 
ferent  geometiy  and  chemical/transport  data.  Nev- 
ertbeless,  the  Markstein  numbers  based  on  Sj  and 
capture  the  general  qualitative  treud  in  terms  of 
(1)  the  leveling'Off  behavior  on  the  rich  side  and  (2) 
the  crossover  poiut  la  the  equivalence  ratio  where 
the  Markstein  number  changes  sign,  wliich  occurs 
la  the  lean  mixture  near  ^  »  0.6  0,7.  On  the  other 

hand  the  Markstein  number  based  on  5^  Is  always 
negative  and  appears  to  be  inadequate  in  capturing 
the  strain  effect  on  the  flame  speed.  This  Issue  has 
also  been  discussed  in  a  previous  study  [25],  Further 
discrepancies  in  5^  caused  by  flow  unsteadiness  are 
discussed  In  the  next  section. 


Response  lo  Oscillalory'  Strain  Rate 

We  performed  a  series  of  calculations  for  pre> 
mixed  flames  with  different  mixture  conditions  sub* 
jected  to  an  osclUatoiy  velocity  field  given  by  equa* 
tion  9.  Table  1  shosvs  the  parametric  values  for  the 
various  test  cases.  Each  case  is  computed  for  various 
frequencies,  mostly  ranging  from  1  to  1000  Hi.  In 
response  to  the  velocity  oscillation,  both  strain  rate 
ajid  flame  speed  fluctuate  In  time,  Afrer  a  number 
of  cycles,  the  solution  approaches  a  limit  cycle  iu  a 
pha^  space  of  flame  spe^  versus  strain  rate. 

Figure  3a  shows  the  result  for  the  <i>  —  0.4  case, 
representing  a  typical  flame  speed  response.  As  the 
frequency  of  oscillation  increases,  responses  of  both 
Sc  and  Kji  are  attenuated,  resulting  In  a  reduced  size 
of  the  limit  cycle  in  Fig.  3a.  Moreover,  the  increased 
phase  delay  of  the  flame  response  to  strain  rate 
causes  tiltuig  of  the  response  cuin«,  resulting  in  a 
smaller  Markstein  number  (the  slope  of  the  ellipse). 
The  result  of  the  flame  response  attenuation  is  con* 
sistent  with  previous  studies  of  unsteady  flames 
[26,27].  Although  not  shown  here,  the  response  of 
Sj  to  unsteady  strain  is  similar  to  either  that  of  S^. 

The  response  of  Su  plotted  in  Fig.  3b,  however, 
shows  drastically  different  results  from  that  of  5^  or 
Sj.  It  appears  as  though  the  flame  speed  oscillation 
is  amplified  with  increasiug  frequency.  Such  a  sub* 
stantlal  difference  can  be  examined  by  comparing 
the  velocity  fields  (or  the  low  ajid  lii^i  frequency 
cases.  Fig.  4  shows  a  comparison  of  the  axial  velocity 
profiles  for  a  quasi>steady  -  1  Hz)  and  a  hig^> 
frequency  response  if  »  500  Hz),  At/  *  1  Hz,  the 
flar^  respomls  to  the  velocity  fluctuation  lu  a  quasi* 
steady  matuier,  such  that  at  any  given  time  the  sc^ 
lution  field  approaches  the  steady  solution  for  the 
instajitaneous  nodle  exit  velocity.  As  a  result,  there 
is  a  larger  lateral  displacement  of  the  flame  location 
and  a  smaller  difference  In  the  local  peak  velocity  as 
u(z  *  0)  varies  from  100  to  200  cm/s.  For  the  hi^- 
frequency  case,  however,  the  reaction  zone  does  uot 
readily  respond  to  the  unsteady  convective  flow. 
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Kwiovitz  nufrber,  Ka 


(a) 


0  1  01S  OJ  09  0» 


KariovC  numbar.  Ke 

(b) 

Fig.  3.  Response  of  (a)  S«  aod  {b)  to  tbe  Kartovitz 
number  Xe  for  various  frequencies  of  osciHaGon;  ^  0.4. 

leading  to  a  smaller  change  in  fiame  location  and  a 
larger  fluctuation  In  5^. 

The  aforementioned  observation  implies  tliat  at 
liigh  frequencies,  tlie  transient  effect  is  balanced 
oiily  by  the  outer  cousectlve  transport,  vdtile  the  in* 
ner  reactive  zone  with  a  shorter  characteristic  time 
scale  remains  largely  unaffected  by  the  unsteady 
oscillation.  In  other  words,  there  is  a  segregation  in 
the  balance  of  the  terms  throu^out  the  flame,  and 
thus  tlie  (ar'Upstream  velocity  cannot  properly  rep¬ 
resent  the  reaction  zone  behavior.  Consistent  results 
were  also  found  in  the  recent  study  of  unsteady 
counterflow  methane/air  premixed  flames  [2&].  In* 
ddentally,  this  explanation  may  account  for  the  ex¬ 
perimental  observation  by  Hlrasawa  et  al  [22], 
where  the  stretch  and  the  flame  speed  fluctuation 
Increase  with  frequency.  In  the  follc^ng  the  results 
are  presented  in  terms  of  5^.  In  favor  of  Its  unique¬ 
ness  in  definition  and  its  physical  relevance  to  the 
overall  burning  rate. 


Axial  aUtanea  (cm| 


Fig.  4.  Companson  of  the  axial  velocity  profiles  for  tbe 
quasi-steady  response  at/  «  1  Hz  (solid  lines)  and  for  tbe 
high  hequeocy  response  at/  «  500  Kz  (dotted  lines],  for 
—  0.4.  Id  each  case,  two  instantaneous  veloaty  flelds  at 
u(x  e  0)  B  100  and  200  cm/s  are  plotted. 


As  derived  from  linearized  theory  [11,12].  the  un¬ 
steady  flame  speed  response  to  strain  rate  can  be 
represented  by  tlie  Marlcsteui  transfer  function,  de¬ 
fined  as 


(U) 


where  the  subscripts  mui  ajid  max  denote  the  min¬ 
imum  and  maximum  values  througli  the  oscillation. 
It  is  expected  tliat  M,.  is  a  fujiction  of  frequency  and 
Lewis  number,  which  depends  on  the  mixture  com¬ 
position  and  the  amoujit  of  dilution.  Tlie  frequency 
is  normalized  by  the  Inverse  of  the  characteristic 
flame  time.  For  flames  wltli  detailed  chemistry  with 
a  number  of  species,  the  proper  defkution  of  the 
characteristic  dame  time  is  not  clear.  One  reasonable 
choice  is  tlie  strain  rate  for  the  initial  stea<^  flame 
for  each  case,  defined  as  Kq,  as  it  properly  represents 
the  characteristic  flow  time  for  the  reference  stea^ 
condition.  The  reference  strain  rate  values  are 
shown  in  Table  1. 

Figure  5  shows  the  normalized  Marksteln  transfer 
function  versus  the  normalized  frequency,  2t^/KQ. 
The  transfer  function  for  each  case  is  normalized  by 
its  quasi-stea<^  value.  As  tlie  frequency  increases, 
most  of  tlie  cuAes  fall  off  at  the  normalized  fre¬ 
quency  of  order  unity.  This  result  implies  that  the 
flame  speed  response  for  the  most  part,  becomes  less 
sensitive  to  a  flow  with  a  busier  transient,  consistent 
with  the  results  from  simplified  theory  [LL,L2]  and 
direct  numerical  simulation  [I3,L4], 

A  unique  response  is  noted,  however,  for  a  range 
of  fuel-rich  conditions,  where  the  transfer  function 
exhibits  a  peak  at  a  moderate  frequency  range  before 
it  eventual^  decays  to  zero.  Such  is  tbe  case  for 
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Fig.  5.  Normalized  MarksCem  transfer  fusctioo  Me 
based  od  as  a  fuocQoo  of  Dormalued  frequeocy, 
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Fig.  6  Response  of  to  (be  Karlovitz  number  ¥Ji  for 
various  frequencies  ofoscillatioD;  (a)  ^  b  3.0^  (b)  ^  *  6.5. 


Fig.  7  Normalized  Marksteio  0‘ansfer  function.  Me 
based  ou  as  a  function  of  normalized  frequency, 
for  ^  B  3.0.  Solid  symbol  denotes  the  undiluted  case,  Og/ 
(O2  +  Ng)  B  0  21,  and  open  symbol  denotes  the  diluted 
case,  Oj/(02  +  Nj  *  0,114. 


•t>  B  3,0,  and  to  a  lesser  e:deQt,  for  ^  ^  5.0,  as 
showa  in  Fig.  5.  To  illustrate  the  cause  for  this  dis> 
tbict  behavior,  Fig,  6a  and  b  show  a  comparison  of 
the  phase-space  response  for  4>  s  3,0  and  6,5.  For 
•t>  —  6.5,  similar  to  the  ^  s  0.4  case  shown  in  Fig. 
3.  the  ellipse  bits  toward  a  direction  such  that  the 
flame  response  attenuates  as  the  frequency  in- 
creases.  On  the  other  band,  for  ^  —  3.0  case,  the 
ellipses  appear  to  be  tilting  up  toward  the  vertical 
direction  for  frequencies  up  to  /  »  2000  Hz.  For 
even  hig^r  frequencies,  the  slope  starts  to  decrease, 
for  example,  at/  »  4000  Hz,  Tliis  su^ests  a  resc^ 
nance  behavior  and  Is  similar  to  the  nonmonotonic 
burning  rate  response  observed  in  a  previous  luiear 
analysis,  as  discussed  by  Im  et  al,  [27],  Similar  be¬ 
havior  was  also  observ^  in  a  recent  study  with  lin¬ 
earized  theory  [12]  where  tlie  peak  in  the  bansfer 
function  was  observed  for  flames  where  the  deficient 
species  has  a  larger  Lewis  number.  The  present  re¬ 
sults  are  distuict  in  that  such  a  peak  response  dis¬ 
appears  again  as  the  equivalence  ratio  further  lu- 
creases,  that  is,  for  even  larger  Lewis  numbers.  On 
the  otlier  bajid,  it  should  be  noted  tliat  an  Increase 
in  the  equivalence  ratio  not  oidy  Increases  the  effec¬ 
tive  Lewis  number  of  the  mixture,  but  also  modifies 
the  flame  temperature  ajid  the  heat  release  param¬ 
eter,  y  *  (/>u  ~  A))''Pu>  which  is  ajiotber  Importajit 
parameter  of  the  problem  in  goveruing  the  Mark- 
stein  transfer  function  [  12] .  Considering  that  tlie  res¬ 
onance  behavior  Is  absent  in  the  earlier  stu<W  with  a 
constant-density  assumption  [11],  such  an  adi^tioual 
effect  may  be  responsible  for  the  nonmonotonic 
trend  with  respect  to  the  equivalence  ratio. 

To  further  assess  the  dilution  effect.  Fig.  7  shows 
a  comparison  of  the  Marksteln  transfer  function  for 
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undiluted  :2nd  diluted  cases  (or  4>  ^  3.0,  The  uadi' 
luted  case  is  the  result  from  Fig.  5  wlille  the  diluted 
case  corresponds  to  air  diluted  with  nitrogen, 

(Og  +  N2)  *  0.114  by  wJume,  such  tliat  the  flame 
speed  and  temperature  are  substajitially  reduced 
(see  Table  1).  The  resonance  behavior  observed  for 
the  ujidiluted  case  is  substanbally  suppressed  with 
dilution,  although  some  nonmonotonic  response  re^ 
mains  at  similar  values  of  the  nonnali2ed  frequency. 
The  results  suggests  that  both  mixture  Lewis  num* 
ber  and  heat  release  parameter  contribute  to  the 
nonmonotonic  flame  response  to  unsteady  strain 
rate. 

Incidentally,  the  fact  the  the  curves  shown  in  Figs. 
5  and  7  are  not  perfectly  smooth  is  probab^  related 
to  the  nonlinear  behavior  of  the  name  speed  re^ 
sponse  to  strain  rate.  For  some  cases  under  study, 
the  range  of  the  strain  rate  fluctuation  is  large 
enough  to  cause  a  substantial  “bending”  of  the  flame 
speed  versus  strain  rate  curve  in  the  steady  reference 
case  {such  as  the  10  Il2  case  in  Fig.  6b).  Therefore, 
the  decrease  in  the  Markstein  transfer  function  with 
frequency  shown  in  Figs.  5  and  7.  which  is  normal' 
i2ed  by  the  steady  reference  case,  may  not  behave 
smoociily  depending  on  the  mean  value  ajid  ampli' 
hide  of  the  imposed  strain  rate  fluctuation.  Never' 
theless.  the  quahtative  results  of  the  paper  remain 
unchanged. 


Conclusions 

The  response  of  hydrogen/air  premixed  flames  to 
osclllatoiy  strain  rates  was  studied  in  an  opposed 
nozsJe  conflguration.  It  svas  found  that,  while  expe' 
dient  for  experimental  measurements,  the  upstream 
local  peak  value  of  the  axial  velocity.  Su,  is  inadequate 
In  representing  the  burning  velocity  characteristics 
for  highly  transient  flames  due  to  the  inability  of  the 
reaction  layer  to  respond  to  the  rapid  transients  in 
the  convective  transport  2one.  Tlie  displacement 
speed,  5^.  is  useful  in  describing  the  flame  as  a  front, 
but  the  results  depend  on  the  specific  choice  of  the 
Isolevel  surface.  The  consumption  speed,  5^  does 
not  suffer  from  the  ambiguity  of  loca^n  and  Is  thus 
maioty  examined  as  a  measure  of  the  integrated 
burmng  velocity. 

The  Markstein  transfer  function  versus  the  fre' 
quency  of  oscillation  was  also  investigated  for 
mixtures  with  various  equivalence  ratios.  The  gem 
eral  trend  is  that  the  flame  speed  response  to  strain 
rate  fluctuation  attenuates  as  the  frequency  exceeds 
the  Inverse  of  the  characteristic  flame  time,  consiS' 
tent  with  many  previous  studies  on  ujisteady  diffu' 
sion/premixed  flames.  Unique  behavior  was  ob^ 
servt^  however,  for  some  rich  mixtures,  where  the 
trajisfer  function  exhibits  a  peak  response  for  mod' 
erate  frequencies.  This  distinct  behavior  disappears 
agaui  as  the  equivalence  ratio  is  further  increased. 


possibly  due  to  the  effect  of  the  reduced  heat'release 
parameter.  Further  comparison  of  tlie  Markstein 
transfer  function  to  the  case  with  dilution  revealed 
that  the  resonance  response  is  suppressed  by  dilm 
don, 
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